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ABSTRACT 

Recent observations of the central compact object in the Kesteven 79 supernova remnant show that 
this neutron star (NS) has a weak dipole magnetic field (a fewx 10^'^G) but an anomalously large 64%) 
pulse fraction in its surface X-ray emission. We explore the idea that a substantial sub-surface magnetic 
field exists in the NS crust, which produces diffuse hot spots on the stellar surface due to anisotropic 
heat conduction, and gives rise to the observed X-ray pulsation. We develop a general-purpose method, 
termed "Temperature Template with Full Transport" (TTFT) , that computes the synthetic pulse profile 
of surface X-ray emission from NSs with arbitrary magnetic field and surface temperature distributions, 
taking into account magnetic atmosphere opacities, beam pattern, vacuum polarization, and gravitational 
light bending. We show that a crustal toroidal magnetic field of order a fewxlO^^'G or higher, varying 
smoothly across the crust, can produce sufficiently distinct surface hot spots to generate the observed 
pulse fraction in the Kes 79 NS. This result suggests that substantial sub-surface magnetic fields, much 
stronger than the "visible" dipole fields, may be buried in the crusts of some young NSs, and such hidden 
magnetic fields can play an important role in their observational manifestations. The general TTFT tool 
we have developed can also be used for studying radiation from other magnetic NSs. 

Subject headings: stars: neutron - magnetic fields - pulsars: individual (PSR J1852-I-0040) - X-rays: 
stars - radiative transfer 



1. INTRODUCTION 

Magnetic field is perhaps the single most important 
quantity that determines the observational manifestations 
of isolated neutron stars (NSs; see, e.g., Turolla 2009; 
Kaspi 2010). Most of our knowledge of NS magnetic fields 
is based on the measurement of the rotation period P and 
period derivative P, which, assuming that NS spin-down 
is due to angular momentum loss through magnetic dipole 
radiation or magnetic wind, yields estimates for the dipo- 
lar magnetic field strength of many NSs. Several NSs are 
observed to have spectral lines which, when interpreted as 
electron/ion cyclotron lines or atomic lines, yield a value 
for the surface magnetic field strength (e.g.. Ho & Lai 
2004; van Kerkwijk & Kaplan 2007; Suleimanov et al. 
2010). Neither of these measurements sheds any light on 
the internal magnetic field configuration. Several recent 
observations suggest that the internal NS magnetic field 
must be more complex than a simple dipole. For exam- 
ple, the pulse profiles of X-ray Dim Isolated Neutron Stars 
may be explained by the addition of higher order multi- 
poles (Zane & Turolla 2006). Also, magnetars are thought 
to require large toroidal fields to fuel their bursting activ- 
ity (e.g., Thompson & Duncan 2001). The most striking 
example of a NS undergoing magnetar-like bursts is SGR 
0418+5729, whose observed dipole field (^ 8 x lO^^G) is 
much too low compared with typical magnetars, and there- 
fore a significant hidden field must be present to power it 
(Rea et al. 2010; Turolla et al. 2011). 

The central compact object (CCO) PSR J1852-I-0040 in 
the Kesteven 79 supernova remnant is another source with 
unusual properties (Halpern & Gotthelf 2010). Specifi- 
cally, while its dipole surface magnetic field is inferred via 
P-P measurement (P = 0.105 s, P = 8.7 x 10"^* s s^^) 
to be only ^ 3.1 x 10"'^° G, its pulse profile presents with 
an unusually high pulse fraction (64 ± 2%), which implies 



large temperature anisotropy (concentrated hot spots) on 
its surface. Such hot spots are unattainable by anisotropic 
heat conduction at the observed dipole field strength. 
Also, since the spindown luminosity is much smaller than 
the X-ray luminosity, they cannot be produced by heat- 
ing from magnetospheric activities. In this paper we ex- 
plore a more complex internal field configuration, adding 
a large crustal toroidal field aligned about the magnetic 
dipole axis to suppress heat conduction everywhere ex- 
cept toward the magnetic poles of the NS. While the idea 
of using crustal field to suppress heat conduction is not 
new (e.g., Geppert et al. 2004, 2006), our goal is to use 
the observed X-ray light curve to quantitatively constrain 
the magnitude and shape of this hidden toroidal field. 

To quantitatively model the X-ray light curve of a 
magnetic NS, it is necessary to add up emissions from 
all surface patches of the star, taking account of beam- 
ing/ anisotropy due to magnetic fields and light bending 
due to general relativity (Zane & Turolla 2006). While 
this is conceptually straightforward, it is not easy to imple- 
ment in practice, as atmosphere models for many surface 
patches (each with different Toff and B) are needed. There- 
fore, in Section 2 we develop a practical method, termed 
"Temperature Template with Full Transport" (TTFT), 
for computing radiation intensity from different surface 
patches. This method captures the essential features of ra- 
diation from magnetic atmospheres, but does not require 
full self-consistent atmosphere modeling for each surface 
patch. The TTFT method allows us to compute NS light 
curves with arbitrary surface magnetic field and tempera- 
ture distributions. In Section 3, we study various crustal 
magnetic field configurations, present the resulting surface 
effective temperature distributions, and apply the method 
of Section 2 to compute the pulse profiles and reproduce 
the observed high pulse fraction of the CCO in Kes 79. We 
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discuss the results in Section 4, and conclude with Section 
5. 

2. LIGHT CURVE CALCULATION AND THE TTFT METHOD 

In this section we outline our method for computing 
X-ray light curves for NSs with arbitrary surface magnetic 
field and temperature distributions. Although our applica- 
tion focuses on the CCO in Kes 79 (Section 3), the method 
presented here is general and can be applied to other types 
of magnetic NSs with atmospheres. □ 

2.1. Light Curve Calculation 

The observed specific radiation flux dF° due to an in- 
finitesimal patch (area dA) on the NS surface is given by 



dF° = lyn" = ^/^(k) cos a 



dcosa 



dfi 



dA, (1) 



where 1° is the observed specific intensity; dil° is the ob- 
served angular size; g is the gravitational redshift factor 
given by g^ = {l-2GM/Rc^)- M, R, and D are the 
mass, radius, and distance to the NS; is the local spe- 
cific intensity (as emitted at the surface patch of the NS); 
k is the direction of photon propagation at the surface; 
a = cos~^(k ■ f) is the angle between k and the local ra- 
dial vector; and fj, is the cosine of the polar angle (with 
respect to the observation axis) of the patch. Due to grav- 
itational light bending, cos a is a non-trivial function of /x. 
Following Beloborodov (2002), we make an approximation 
to this function, and calculate the observed specific flux as 



F° 



9 

Z?2 



/^(k) [l-g\l-^i)\dA, 
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where the integration is over the entire observable area (in- 
cluding light bending) of the NS surface. Thus, given the 
knowledge of local radiative intensity [/^(k)] of all surface 
patches, combined with the photon trajectory at differ- 
ent rotational phase (time) [/i(t) <^ Q!(t)], the resultant 
observed specific flux can be calculated. 

2.2. Radiative Transfer in Magnetic Neutron Star 
Atmospheres 

The full computation of /^(k) requires construction of 
atmospheric models in radiative equilibrium for various 
surface magnetic fields B (magnitude and direction), effec- 
tive temperatures Toff, and compositions (e.g., Shibanov et 
al. 1992; Pavlov et al. 1995; Zane et al. 2001; Ho & Lai 
2001, 2003, 2004; Potekhin et al. 2004; van Adelsberg & 
Lai 2006; Suleimanov et al. 2010). This is currently not 
practical for exploration of a wide range of model param- 
eters, since quantitative calculation of the light curve de- 
mands adding up many different surface patches (rather 
than one or two hot spots; cf. Ho 2007). Therefore we 
develop an approximate method (TTFT) (see Section 2.3) 
that allows us to obtain /^(k) efficiently. Here, we describe 
our full radiative transfer calculation in NS atmospheres. 

In the atmospheres of magnetized NSs, photon polar- 
ization plays an important role in determining /^(k). For 
example, at photon frequencies smaller than the electron 
cyclotron frequency, the mean free path for ordinary (O) 

^ We do not consider the extreme case of high magnetic field at low 
Lai & Salpeter 1997; Medin & Lai 2007). 



mode photons is much smaller than that of the extraordi- 
nary (X) mode photons, and therefore, the X-mode pho- 
tons that contribute to the emergent flux are typically pro- 
duced in the deeper, hotter layers of the atmosphere. At 
higher field strengths (> 8 x lO^^G), vacuum polarization 
will also have a strong effect on the radiation spectrum 
(Lai & Ho 2002, 2003). While rigorous treatment of this 
effect is available (van Adelsberg & Lai 2006), we describe 
an approximate method below. 

For a photon of a given energy E propagating along k 
in the NS atmosphere, vacuum resonance occurs at 



pv = o.96y. 



'^ElBlJs^ g cm- 



(3) 



where is the electron fraction, Ei — E/{1 keV), B14 = 
i?/(10^^ G), and Jb is a slowly varying function of B, 
of order unity. At resonance, the X and O polarization 
modes become degenerate, and therefore the photon may 
be converted from one mode to another when it traverses 
the resonance. We define two more effective polarizations 
modes: -I- and -. These modes act as X (O) polarization 
modes before the resonance and as O (X) modes after. Lai 
& Ho (2002) showed that photons satisfying the adiabatic 
condition, E > i?ad, act as -I- and ~ modes. Here 
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keV, 
(4) 

Ebz = OMYeBii keV is the ion 
= \ds/dhip\ is the density scale 
height along the ray. More generally, a photon with energy 
E ^ -Ead will have a nonadiabatic 'jump' probability of 



where 9kB = cos ^ k • B, 
cyclotron energy, and Hp 



Pi, 



exp 



-^iE/E,,f 
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If the 'jump' probability is low, most outward diffusing 
photons will evolve through the resonance adiabatically 
as -|- or - modes, undergoing mode conversion. On the 
other hand, if the 'jump' probability is high, X-mode (O- 
mode) photons will remain X-mode (0-mode) photons af- 
ter traversing the resonance. We therefore elect to treat 
the X, O, and - photon modes as distinct carriers whose 
relative abundances are determined by Pjump: 

i: = Pjumpil^ + I?) + (1 - Pjump)(/+ + /-). (6) 

Note that, for the application to the Kes 79 CCO (Section 
3), at fields B ~ lO^'^ - 10" G and energies E ^l- 5keV, 
Pv is smaller than the photosphere densities, and therefore 
vacuum polarization does not play an important role. We 
therefore treat all photons as "jump" photons (Pjump = 1)) 
and do not compute intensities for the -I- and - modes. 

In general, to compute the individual specific intensities 
for each carrier, we integrate the standard radiative trans- 
fer equation, dll/drl = B^ — II. However, for the present 
purposes it is sufficient to use the Eddington-Barbier rela- 
tion: 



Ilik) « ^B, [Tirlik) = 2/3)] 



(7) 

temperature, for which the NS surface may be in condensed form (e.g., 
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Fig. 1. — Top panel: NS atmosphere temperature profiles vs 
Thomson depth from several self-consistent atmosphere models. 
Solid hne - B = 6.2 X lO^" G, T^g = 4 X lO'' K; dashed line 
- B = 6.2 X IQl" G, Tcff = 1 X 10^ K; dotted line - B = 0, 
Teff = 5 X 10^ K; dashed-dotted line - B = 0, T^b = 1-6 X 10*^ K; 
dashed-triple-dotted line - S = 2 x lO^^ Q, T^s = 5x 10^ K. Bottom 
panel: The same profiles rescaled such that Tgfi = 3.8 x 10® K at 
a normal magnetic field with strength B = 8.8 x 10^° G using the 
TTFT method. In addition, we include the Thomson grey opacity 
temperature profile ~ (r + 2/3) (long dashed), similarly rescaled. 

That is, the emergent intensity (toward k) due to carrier 
i is equal to the Planck function cvahiatcd at the depth 
at which the optical depth for that carrier (along k, 
measured from the outside in) is equal to 2/3. The fac- 
tor of 1/2 arises because the photon number distribution 
is split among two carriers at all times. Thus, to calcu- 
late the total emergent specific intensity, only a few pieces 
of information are required. These are: the direction and 
magnitude of the magnetic field, the opacities for each car- 
rier, and the atmospheric temperature and density profiles 
[T(r) and p{t), where r = ktV (with kt = 0.4 cm^g~^) 
is the (zero-field) Thomson depth, and y is the column 
density]. 

We use free- free absorption opacities for each carrier, 
inchiding the contribution of electron-ion Coulomb col- 
lisions, as described by Lai & Ho (2003). This opacity 
depends strongly on k ■ B, as well as on density and tem- 
perature. 

In addition, at the magnetic field strengths we are con- 
sidering, the opacity is further complicated by oscillations 
in the magnetic Gaunt factor (Potekhin 2010). However, 
the energies we consider have values of ~ 1.5 — 5 times the 
electron cyclotron energy and therefore these oscillations 
are not very strong. Furthermore, while the oscillations 
may affect the specific flux emitted at certain frequencies, 
we do not expect them to affect pulse shape. We there- 
fore do not include this complication in our calculation, 
instead choosing to use the non-magnetic Gaunt factor, 
which captures the trend about which the full magnetic 
Gaunt factor oscillates. 

The temperature profile T(r) and density profile p(t) 
are related by the equation of hydrostatic equilibrium. For 
all practical purposes, an ideal, non-degenerate gas equa- 
tion of state can be used with better than a few percent 
accuracy (e.g., Ho & Lai 2001). 
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Fig. 2. — Sample full-star light curves of NS surface emission at 
photon energies around 1.5 keV calculated using the TTFT method. 
Two rotation phases are shown. These light curves are computed 
for a specific NS surface temperature and magnetic field distribu- 
tions and viewing geometry (see Section 3). The different curves 
correspond to the results obtained using different temperature pro- 
file templates as depicted in Figure 1. The highest pulse fraction 
obtained is 62%, and the lowest 53%. 

2.3. TTFT: A Method for Temperature Profile 

Computation 

As an all-purpose tool, the computational procedure 
described in Sections 2.1-2.2 requires the knowledge of 
the atmospheric temperature profile T{t). In theory, for 
a given combination of effective temperature, magnetic 
field strength and direction, T(r) is obtained by self- 
consistent atmosphere modeling to achieve radiative equi- 
librium. However, this modeling must be done for every 
single patch of the NS surface, and is not practical to im- 
plement. Here we propose the "Temperature Template 
with Full Transport" (TTFT) method for constructing 
temperature profiles. While this method is not accurate 
for predicting the observed synthetic radiation spectrum, 
we show that it is reliable insofar as light curve shape (and 
therefore pulse fraction) is concerned. 

The TTFT method consists of the following steps: (i) 
We write the temperature profile in the form T(t) = 
ToF(t), characterized by the dimensionless template func- 
tion F{t) and the adjustable constant Tq. The template 
function may be obtained from a specific NS atmosphere 
model, (ii) For each T(t), we carry out full radiative trans- 
fer calculation (Section 2.2) to obtain the total radiation 
flux from a single surface patch, 

T= Jdujd'^]<il^{k)\i-r, (8) 

where the integration / rf^k covers the outgoing direction 
relative to the surface normal vector r. Since T — (yT^ff, 
we then obtain a one-to-one mapping between Tq and 
Teff for the given template function F{t). (iii) We then 
adopt the ansatz that for a given surface patch with T^s, 
the atmosphere profile is T(r;Teff) = To(Teff)F(r), and 
the emergent radiation intensity can be computed via the 
method of Section 2.2. 

Figure 1 demonstrates the reasoning for the validity of 
the TTFT method. The top panel presents several self- 
consistent temperature profiles at different effective tem- 
peratures and field strengths provided by Wynn Ho {B = 
and B = 2 X 10^^ G models) and by Valery Suleimanov 
{B = 6.2 X 10^° G models). We note that all the profiles 
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are very similar in shape, especially for r <; 0.1, and dif- 
fer mostly in overall normalization; this is true both with 
respect to magnetic field strength variation and effective 
temperature variation. After normalization to the same 
Teff = 3.8 X 10^ K and B = 8.8 x 10^° G, ah the temper- 
ature profiles become similar at t 0.1. Therefore, when 
computing the intensity of photons that emerge from suf- 
ficient depth in the atmosphere, the same single tempera- 
ture template F{t) may be used for the entirety of the NS 
surface. 

Note that it is not surprising that all the temperature 
profiles are essentially identical at t > 1; indeed, in the 
diffusion-limited regime that is to be expected. The less 
obvious and most crucial point is that the profiles remain 
very similar for 0.1 < t < 1, where they also differ signifi- 
cantly from a grey opacity profile. In our computations of 
synthetic light curves (see Figure 2), we see that only the 
photons traveling nearly radially out from the star emerge 
from slightly deeper than r ~ 1, and thus most of the pho- 
tons contributing to the light curve emerge from smaller 
Thomson depths. Therefore, a grey opacity temperature 
profile could not possibly produce the same large synthetic 
pulse fraction, and would not yield a realistic representa- 
tion of the emergent flux; we must use a self-consistent 
template accounting for depth-dependent opacity. 

Figure 2 shows sample X-ray light curves (normalized 
to a mean of unity) of magnetic NS surface emission, 
computed using our TTFT method with various template 
functions F{t) depicted in Figure 1. These light curves 
are for a specific NS surface model (effective temperature 
and magnetic field distributions) appropriate for the CCO 
in Kes 79 (see Section 3 for details). We see that the 
shapes of these light curves are very similar, with the pulse 
fraction ranging from 53% to 62%, despite the fact that 
the template functions are obtained from vastly different 
atmosphere models. Note that in renormalizing the tem- 
perature profile (to obtain the Tq-Tcb mapping), we make 
the approximation that the magnetic field is normal to 
the stellar surface, and we keep its magnitude always set 
to the field strength at the magnetic poles. This speeds 
up the construction of the temperature profile (i.e., the 
2o-Toff mapping for a given template function) since the 
integration in Equation (|8]) can be done more efficiently, 
and need not be repeated for every surface patch. We have 
found that this approximation introduces rather small er- 
ror in the resulting Tq-Tcs relation; this is understandable 
because T^s is determined by summing over the emitted 
spectrum in all directions, and the direction of the mag- 
netic field is averaged over, while the radiative transfer is 
not too sensitive to factor-of-two changes in field strength. 
Nevertheless, we note that when computing /^(k) for the 
light curve, the "correct" magnetic strength and direction 
must be adopted in order to accurately account for the 
magnetic beaming effect. 



3. CONSTRAINING THE CRUSTAL MAGNETIC FIELD OF 
THE NEUTRON STAR IN KES 79 

We now use the method described in Section 2 to model 
the X-ray light curve of the NS in Kes 79, thereby con- 
straining the magnetic field configuration in the NS crust. 




Fig. 3. — Effective surface temperature distribution T^f{{0) 
(heavy lines) for different crustal magnetic field configurations. The 
toroidal field is characterized by Equations JOll-l llOI I. Solid line — 
Bo = 2 X IQl* G, do = 40°; dashed line - Bo = 2 X 10^** G, 
Oo = 10° ; dot-dashed line - Bo = 2 X lO^^ G, 6»o = 40° . The light 
lines show the corresponding function f{d) at different 9o- 



An internal field of strength similar to the inferred dipole 
field of the Kes 79 NS would lead to a very uniform tem- 
perature distribution on the NS surface (e.g., Schaaf 1990; 
Potekhin et al. 2003). In order to produce sufficiently con- 
centrated hot spots, we introduce a toroidal (with respect 
to the dipole axis) component to the magnetic field, such 
that it is sufficiently large everywhere except in a concen- 
trated region around the magnetic axis. We adopt the 
following form for the toroidal field: 



B^{r,e) = B^f{e)g{r) 



(9) 



where Bq is an adjustable maximum field strength and 9 
is the magnetic colatitude. 

Because the NS atmosphere is likely too tenuous to be 
able to sustain a toroidal field, the toroidal field should 
be contained entirely in the crust, such that the poloidal 
current producing it is also contained in the crust. Thus, 
a realistic form of g{r) should vanish at least at the NS 
surface, and possibly at the crust-core interface as well. 
However, for simplicity, in the present calculation we treat 
g{r) as constant. We choose the following functional form 
for /: 



f{0) = 



/o(sin0/sin0o)' 



1 



(sine/ sin 00)2' 



(10) 



where 9o is an adjustable parameter defining how sharply 
the toroidal field rises away from the magnetic poles. The 
normalization constant /o is specified so that max[/(0)] = 
1. Obviously, this functional form for f{6) is representa- 
tive; other forms are possible. However, the detailed shape 
of f{9) cannot be realistically constrained and is not im- 
portant for our purpose. Our goal is to examine what 
ranges of values of Bo and are required in order to pro- 
duce sufficiently distinct surface hot spots to generate the 
observed large pulse fraction in Kes 79 NS. Note that we 
are essentially exploring a parameter space orthogonal to 
the one considered by Geppert et al. (2006), who fixed 
f{9) to be a simple sinusoid and explored the dependence 
of surface temperature and pulse shape on g(r). 
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Fig. 4. — Top panel: Full-star synthetic light curves at photon 
energies around 1.5 keV for different crustal magnetic fields. The 
different lines correspond to the different effective temperature dis- 
tributions Tofi-(6) depicted in Figure 3. The pulse fractions are 58% 
(solid line), 55% (dashed line) and 38% (dot-dashed line). Bottom 
panel: Light curves for the Sq = 2 x 10^"' G, 6o = 40° model at 
different photon energies: E ~ 1.5 keV (dotted line, pulse fra<;tion 
58%), B ~ 3 keV (solid line, PF 58%) and ~ 5 keV (dashed line, 
PF 63%). 

In addition to the crustal toroidal magnetic field, we also 
include a star-centered dipole field. The canonical surface 
dipolc field strength (at the magnetic poles) as inferred 
from the P-P measurement is Bp sin ~ 6.2 x 10^° G 
(for R = 10 km and moment of inertia ~ 10^® g cm^), 
where 9b is the angle between the magnetic axis and the 
spin axis. Of course, pulsar spindown power is not ex- 
actly given by the vacuum dipole radiation formula (Con- 
topoulos & Spitkovsky 2006). For concreteness, we adopt 
Bp ~ 8.8 X 10^° G in our calculations (corresponding to 
9b = 45°), although using a somewhat different value 
(within a factor of two) would only have a small effect 
on our result. 

Given the magnetic field configuration, the effective 
temperature distribution on the NS surface is determined 
by heat conduction in the crust (e.g., Schaaf 1990; Heyl 
& Hernquist 1998; Potekhin & Yakovlev 2001) Since the 
magnetic field is symmetric about the magnetic axis, the 
effective temperature depends only on the polar angle, 
Teff = Teg{9). We use the zero- accretion fitting formula of 
Potekhin et al. (2003) to calculate Teg{9). This formula 
takes as inputs a core temperature and a crustal magnetic 
field of arbitrary strength and direction (assumed constant 
through the radial extent of the crust), and returns the ef- 
fective surface temperature. We fix the core temperature 
at Tc = lO^K, corresponding to T^s = 3.8 x 10^ K at the 
magnetic poles, which matches sufficiently the hot spot 
temperature expected from blackbody fits. 

Figure 3 presents three sample ftmctions Toff (0), two at 
the same value of Bq but different values of and one at a 
lower value of Bq. We also include the corresponding f{9) 
curves to clarify the effect of changing 9o. It is important 
to note that when the toroidal field strength Bq is suffi- 
ciently large, even a very smooth B^ profile can produce 
rather small hot spots. For example, at i?o = 2 x 10^^ G 
and ^0 = 40°, the temperature drops from Tp ~ 4 x 10^ K 
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Fig. 5.— PuU-NS light curves at Bq = 2 x lO" G, = 40°, 
E = 1.5 keV, and different geometries. Solid - /3 = 80°, 0b = 90°. 
Dotted - 13 = 80° , Ob = 85°. Dashed - /3 = 80°, 6»b = 80°. 

at the magnetic pole to Tp/2 at 9 ^ 3°, and to Tp/4 at 
9 ~ 10°. 

When computing the synthetic X-ray light curves, we 
split the effective temperature curve into 7-9 bins, with 
degree-resolution close to the pole, and coarser resolution 
for the rest of the star. For each bin we create the tempera- 
ture profile T{t; Teg) at the; requisite effective temperature 
using the TTFT process described in Section 2. Note that 
since B^ vanishes at the NS surface, it does not contribute 
to the atmospheric physics in any way. We can then com- 
pute the light curves for each of the effective temperature 
distributions depicted in Figure 3. 

In addition to Bq and ^o, the light curves also depend 
on two geometrical parameters: the angle between the line 
of sight and the rotation axis, /3, and the angle between 
the magnetic dipole axis and the rotation axis, 0b- We 
select 13 and in order to maximize the pulse fraction. 
Since, as a rule, emission is strongest for photons emitted 
parallel to the direction of the magnetic field, we should 
select P and 9b to be the same or nearly so, such that at 
some point during the NS's rotation, the photons we see 
are emitted parallel to the magnetic field and the beam- 
ing effect is maximized. We also select the angles such 
that a single-peaked pulse profile is produced (see below 
for double-peaked profiles). After exploring the available 
phase space, we settle on (3 = 9b = 45°. This choice is not 
an overly specific one, however, as there are many other ge- 
ometrical configurations that result in only slightly smaller 
pulse fractions. Our selection of /3 and 9b maximizes the 
pulse fraction, and therefore minimizes how "dark" the 
stellar surface must be away from the hot spots. This se- 
lection is therefore the most conservative estimate of the 
required strength of B^. 

The top panel of Figure 4 presents the whole-star syn- 
thetic light curves at the photon energy E = 1.5 keV, 
using each of the three T'eff(^) functions of Figure 3. The 
So = 2 X 10" G, 6*0 = 40° model (soHd fine) achieves 
the required pulse fraction observed in Kes 79 NS. For 
Bo = 2 X 10" G, 6^0 = 10° model (dashed line), the 
hot spot is too concentrated with a small area, caus- 
ing a decrease in pulse fraction. For Bq = 2 x 10^^ G, 
9o = 40° model (dot-dashed line), the spot has enough 
area, but the rest of the star is simply not dark enough, 
and thus the pulse fraction is small. We therefore con- 
clude that a toroidal magnetic field in the crust of order 
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~ a few X 10^^ G or higher, rising slowly enough from 
zero toward the magnetic equator, such that it allows the 
hot spot to gain enough area, can naturally produce the 
requisite pulse fraction observed in the Kes 79 NS. 

The bottom panel of Figure 4 shows that the pulse pro- 
files at different photon energies are approximately the 
same (with the pulse fraction slightly higher at E = 5 keV 
than at 1.5 keV), in agreement with the observation of the 
Kes 79 NS (Halpern & Gotthelf 2010). Note that while we 
compute the pulse profiles at a single well-defined energy, 
they are nevertheless representative of the pulse profiles 
that would be observed in energy bins centered on the 
quoted energies, with widths of a few tenths of a keV. 

Although our theoretical light curves reproduce the 
large pulse fraction observed in the Kes 79 NS, there re- 
mains appreciable difference between the shape of the light 
curves presented in Figure 4 and the observation. Our 
model light curves have relatively sharp peaks, with the 
system spending nearly half of the rotation period at the 
(observed) flux minimum. This is different from the broad 
flux maximum observed in Kes 79 NS. Note, of course, that 
we have only considered one particular functional form for 
the magnetic field distribution, and the form we consider 
is very smooth. It is not clear that more complex mag- 
netic field configurations would yield the right pulse shape 
while maintaining a high pulse fraction. Another possible 
way to resolve this difference is to consider pulse profiles 
with two peaks per rotation period, one peak for each mag- 
netic hot spot. This is possible when both the line of sight 
and the magnetic dipole axis are nearly perpendicular to 
the rotation axis. Such a light curve would be observation- 
ally indistinguishable from a single-peaked light curve, but 
would imply that the actual rotation period of the NS is 
two times shorter. 

Figure 5 presents several examples of the two-peaked 
light curves for different combinations of the geometric 
parameters 13 and 9b (both around 90°). We see that 
the width of the peaks is much closer to what is observed 
in Kes 79 NS. Obviously, when either (3 ov 9b deviates 
significantly from 90°, the two peaks will have unequal 
heights, which can be measurable. Still, our result in Fig- 
ure 5 shows that the square of geometrical phase space 
with corners at {(3,9b) = (80°, 80°) and (by symmetry) 
(110°, 110°) provides viable two-peaked light curves that 
may be indistinguishable from single-peaked light curves. 

4. DISCUSSION 

4.1. Large Pulse Fraction - How? 

Beloborodov (2002) demonstrated that two infinitesimal 
antipodal hot spots, emitting isotropically, on the surface 
of a 1.4Msun, 10 km NS cannot, in full geometric general- 
ity, produce light curves with pulse fractions higher than 
about 12%. Poutanen & Beloborodov (2006) later gener- 
alized this to anisotropic emergent intensities, producing 
pulse fractions of a few x 10%. Thus, it seems that even 
if one were to find a way to create hot spots on the NS 
surface, it would be difficult produce the 60% pulse frac- 
tion observed in Kes 79 NS. How, then, did we accomplish 
this? 

In fact, diffuse hot spots of finite size (with varying tem- 
peratures), combined with the beaming due to anisotropic 
photon opacities in magnetic fields, are essential for gen- 



erating the large pulse fraction in our model. In our cal- 
culations, the hottest regions with radii of ^ 0.5° around 
each magnetic pole do not produce pulse fractions higher 
than about 40%. Instead, it is the subsequent annuli of 
decreasing effective temperature that create the bulk of 
the pulsed emission. Therefore, to quantitatively model 
the light curves of NS surface emission, it is important 
to go beyond the single-temperature hot spot approxima- 
tion widely used in many previous studies. It is for this 
reason that a reasonably efficient method, such as TTFT 
developed in Section 2, is necessary in order to add up 
atmosphere emissions from many different patches of the 
NS surface to produce synthetic light curves. 

4.2. Crustal Magnetic Field Evolutionary Timescale 

The crustal magnetic field assumed in our model (Sec- 
tion 3) is unlikely to be unique for producing the observed 
light curve of Kes 79 NS. In a way, it is the most conserva- 
tive magnetic configuration, with the minimum number of 
free parameters, and allows for self-consistent calculation 
of the surface temperature distribution and the resulting 
light curve. One can invoke more complex magnetic field 
topology, but the general conclusion of our result, that the 
NS in Kes 79 contains a significant hidden magnetic field 
(of order a fewxlO^'^G or larger) in the bulk region of its 
crust, is likely to be robust. We also note that while the 
"strong toroidal field -|- weak dipole field" configuration 
adopted in our model has generally been thought to be 
unstable in fluid stars (e.g., Braithwaite & Nordlund 2006; 
Reisencgger 2009), this may not, in fact, be the case, as 
Braithwaite (2009) demonstrated that the effect of gravi- 
tational stratification is such that even a relatively weak 
poloidal field can be enough to stabilize the configuration. 
In NSs beyond the fluid stage of evolution, the crust also 
has a stabilizing effect. If we include a poloidal field of sim- 
ilar strength as the toroidal field, our general conclusion 
will be unchanged. 

While it is beyond the scope of this paper to specu- 
late how the magnetic field configuration adopted in our 
model may come about, we can estimate how long this 
magnetic field configuration lasts in its current form to al- 
low for a good chance of being observed. Given the toroidal 
field strength we have found above for Kes 79 10^'^ G), 
its evolution should be dominated by Hall drift, since the 
Ohmic timescale for a young NS is much longer than the 
Hall timescale (see Reisenegger et al. 2005; Pons & Gep- 
pert 2007; Pons et al. 2009): 

'dt 



— = -V X 



(V X B) X B 



(11) 



AirneC 

where rig is the electron number density. The characteris- 
tic Hall timescale is 

tHall = , (12) 

where L and B arc the typical length scale and field 
strength. Using values rig ~ lO'^^cm"'^, L ~ 0.5i?, with 
H the crust thickness ~ 1 km, and B ~ 2 x 10"^^ G, 
we find tHaii ~ 8 kyr. We can further evaluate the 
Hall timescale as a function of location in the crust via 
iHaii ~ lB|/|9B/9i|. We find that the Hall time scale is 
greater than 10 kyr throughout most of the crust, and de- 
creases to a value of 5 kyr at the edges of the crust. Thus, 
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the magnetic field configuration is stable enough to remain 
virtually unchanged up to more than the current age of the 
NS (about 7 kyr). 

5. CONCLUSION 

This paper has two main results: 

(i) We have developed a general-purpose method ( "Tem- 
perature Template with Full Transport", TTFT) to ef- 
ficiently compute the radiation intensity from different 
patches of a NS surface with arbitrary magnetic fields 
and effective temperatures. The TTFT method accounts 
for the effects of realistic polarized photon opacities and 
anisotropic radiative transport in NS atmospheres with- 
out self-consistent atmosphere modeling for every surface 
patch. This method is therefore useful for computing syn- 
thetic light curves of X-ray emission from the whole neu- 
tron star surface. 

(ii) We have used the TTFT method to model the X- 
ray light curve of the young NS in the Kes 79 supernova 
remnant. We find that a crustal magnetic field of order 
a fewxlO^'^G, much stronger than the observed dipole 
field, can produce significant surface temperature inho- 
mogeneity and generate the observed large pulse fraction. 
Small hot spots can be produced even for fields that vary 
smoothly from the magnetic poles to the equator. Al- 
though our specific crustal field model adopts a conserva- 
tive axisymmetric toroidal configuration, we believe that 
our general conclusion about the strength of the hidden 
magnetic field in the star's crust is robust for more gen- 
eral field configurations. 

Although we have developed the TTFT method with the 
application to the Kes 79 NS in mind, the basic method- 
ology may prove useful for wider applications. In many 
published works on surface emission from magnetic NSs, 
the observed radiation is compared to a single-temperature 
NS atmosphere model, while in reality, a synthetic radia- 



tion flux from the whole star is needed. Using the TTFT 
method, the total radiation from different patches of the 
NS may be calculated with temperature templates from 
a small number of self-consistent atmosphere models. We 
have shown in this paper that at least with respect to the 
light curves, the TTFT method can produce accurate re- 
sult. 

Our finding that the Kes 79 NS has a much larger hid- 
den magnetic field in the crust than the dipole field in- 
ferred from the spindown measurement (see Halpern & 
Gotthelf 2010) may have important implications for the 
nature of central compact objects (CCOs) and the evo- 
lution of NS magnetic fields. CCOs constitute a signif- 
icant population of NSs in supernova remnants. Recent 
observations have shown that the 7-10 known CCOs (or 
candidates) have weak dipole magnetic fields (Halpern & 
Gotthelf 2010,2011), so their relationship to other types 
of isolated NSs (radio pulsars, magnetars and dim iso- 
lated NSs) is unclear (van Kerkwijk & Kaplan 2007; Kaspi 
2010). The presence of a large sub-surface magnetic field 
in the Kes 79 NS supports the notion that supernova fall- 
back may have submerged the magnetic field, which may 
later diffuse back to the surface (e.g., Muslimov & Page 
1995; Geppert et al. 1999; Ho 2011). Together, this prop- 
erty of the Kes 79 NS and the recent discovery of SGR 
0418-f 5729 with a weak dipole field (Rea et al. 2010; Tur- 
oUa et al. 2011) suggest that both the hidden sub-surface 
magnetic field and the "visible" dipole field of NSs play an 
important role in their observational manifestations. 

We thank Wynn Ho and Valery Suleimanov for sup- 
plying their neutron star atmosphere models, which are 
essential for calibrating our method. We also thank Wynn 
Ho, Alexander Potekhin and Marten van Kerkwijk for 
useful discussion and suggestions. This work was sup- 
ported in part by NSF grant AST-1008245 and NASA 
grant NNX10AP19G. 
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